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ABSTRACT 



The Lyman alpha line is a robust tracer of high redshift galaxies. We present estimates of Lyman alpha emission from a proto- 
galactic halo illuminated by UV background radiation fields with various intensities. For this purpose, we performed cosmological 
hydrodynamics simulations with the adaptive mesh refinement code FLASH, including a detailed network for primordial chemistry, 
comprising the formation of primordial molecules, a multi-level model for the hydrogen atom as well as the photo-ionization and 
photo-dissociation processes in a UV background. We find that the presence of a background radiation field J 2 i excites the emission 
of Lyman alpha photons, increasing the Lyman a luminosity up to two orders of magnitude. For a halo of ~ 10 10 M Q , we find that 
a maximum flux of 5 X 10~ 15 erg cm -2 s _1 is obtained for J 2 \ X f esc =0.1, where f esc is the escape fraction of the ionizing radiation. 
Depending on the environmental conditions, the flux may vary by three orders of magnitude. For J21 x f esc > 0.1 the Lyman alpha 
luminosity decreases as the atomic hydrogen abundance becomes rather small. The fluxes derived here can be probed using Subaru 
and the upcoming James Webb Space Telescope. The emission of Lyman alpha photons is extended and comes from the envelope of 
the halo rather than its core. In the center of the halo, line trapping becomes effective above columns of 10 22 cm -2 and suppresses 
the emission of Lyman alpha. In addition, cooling by primordial molecules may decrease the gas temperature in the central region, 
which further reduces Lyman a emission. In the central core, H2 is photo-dissociated for a background flux of J21 > 1000. For weaker 
radiation fields, i.e. J 2 i < 0.1, H 2 and HD cooling are particularly strong in the center of the halo, leading to gas temperatures as low 
as ~ 100 K. We also performed a parameter study with different escape fractions of ionizing photons and explored the relative role of 
ionizing and dissociating radiation. We find that Lyman alpha emission depends more on the strength of the ionizing background. For 
a constant ionizing background, the Lyman a flux increases at least by an order of magnitude for stronger photodissociation. 
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1. Introduction 



Cosmology: theory - early Universe - Galaxies: formation - Atomic processes - Molecular 



The Lyman alpha line is an important probe of the first 
galaxies as they are potentially strong Lyman alpha emitters. 
Recently, observations of galaxies have been reported at z > 7 
dBouwens et al.ll201ll: iBunker et aDl2010l: IVanzella et ai1l2010l) . 
The upcoming James Webb Space Telescope (JWST) will fur- 
ther enhance our ability to detect galaxies at higher redshift by 
exploring them with higher sensitivity and test our current un- 
derstanding of early structure formation. The detection of the 
Lyman alpha line will be a key feature in the observation of 
galaxies in the early universe. 



pha photons. This scenario has also been supported by theo- 
retical studies and numerical simulation s ( Haiman et al J I200C 
Diikstra et all 120061* lYang et al.l 120061: IDiikstra et al.1 12006a 



Wise et al] 120081: IDiikstra & Loebl 12009b iDekel et al] 1200 



far at intermediate redshifts (Lehnert et al. 20101; 


Steidel et al. 


2000; 


Ouchietal. 2009; Saitoetal. 2008, 2006 


; lYans et al. 


2009; 


Matsuda et al. 20041). The origin and source of these 



blobs is still not completely known. The energy sources of 
these Lyman alpha blobs include photo-ionization by stars o r 
miniquasars (lHaiman & Reesl 2001; J imenez & Haimanll2006l) . 
Another potential driver of these emitters could be the emis- 
sion from cold stream s and accretion flows dNilsson et al.l l2006; 
I Smith & Jarvisl 12007b . The accreted gas releases its binding en- 
ergy in the gravitational potential of the halo, which results 
in the generation of shocks and the emission of Lyman al- 



Goerdt et al] 1201(3; iPawlik et al.1 1201QI: iFauc her-Gig uere et all 
120101; lLatif et afll2011bl) . Some of the observed Lyman alpha 
blobs have been asso ciated with massive star forming galaxies 
dMatsuda et al.ll2006l) . 

In our previous study (lLatif et al] 12011 ah . we estimated the 
emission of Lyman alpha photons from a protogalactic halo 
in the absence of UV radiation. We found that the emission 
of Lyman alpha radiation is extended and originates from ac- 
cretion flows and virial shocks. The situation is however more 
complicated, due to the presence of ionizing and photodisso- 
ciating backgrounds , as well as the formation of molecules 
(I Johnson et al.ll2007l) . These are the effects we explore here. 

The intense Lyman- Werner background flux can photo- 
dissociate H2 molecules in atomic cooling halos with virial 



temperature > 


► 10 z 


1 K (Machaceketal.11200] 


Lb Omukail 120011: 


Johnson et al.l 


2007 


:IO'Shea & Normanll2008l: 


Greif et aDl2008t 


Omukai et al.1 


2008 


; Diikstra et al.1 2008; Shans^taLl 20101). It 



is found that the strength of the UV flux below J21 = 10 3 
onl y delays the collapse and does not quench the formation of 
H? . iSchleicher et ID (1201 Obi) confirmed the previous results and 
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found that for J21 > 10 3 H2 does not form, but the gas tempera- 
ture nevertheless declines down to 5000 K due to cooling from 
higher electronic transitions of hydrogen atom. 

The implications of ionizing radiation for the chemistry 
in the halo are, on the other hand, much more unexplored. 
During the epoch of reionization, background UV radiation gen- 
erat ed by stellar populations will photoionize the intergalactic 
gas ( Gnedin 2000; Schleicher et al. 2008a; Thomas & Zaroubil 
1200a IWise & A bel 2008a; Maselli et^|2009|). The recombi- 



nation following the absorption of stellar photons may fur- 
ther boost the Lyman alpha emission from accretion flows. 
iFurlanetto et al.l ([2005) found that the photoionizing background 
can power the Lyman alpha emissivity and can explain the ob- 
served Lyman alpha blobs. Our main goal here is to study the 
impact of background UV on Lyman alpha emission and to see 
how it influences the spatial emission profile. 

Gas ionized by UV radiation cools faster than it recom- 
bines, preserving an increased ionization degree. The free elec- 
trons act^_as_cjrt^riyjrts i 

gendSaslaw & Zipovl [l967l: iGalli & Palial 119981: lYoshida et al.l 
120071) . In the absence of dust, H2 formation proceeds via the gas 
phase, with the dominant formation channel being due to H" : 

H + e-^lT+y. (1) 

This then leads to H2 formation via the reaction 



H + H" -> H 2 + e" 



(2) 



Another prominent coolant in primordial gas is HD. Its dom- 
inant formation path is via collisions between H2 molecules and 
ionized deuterium: 



H 2 + D + -> HD + H + 



(3) 



While H2 molecules can cool the gas to a few hundred Kelvin, 
HD can decrease the temperature even further if it is abundant 
enough. While the H2 molecule is perfectly symmetric and thus 
has no permanent dipole moment, the deuterium breaks the sym- 
metry of the molecule, thus creating a small but non-negligible 
dipole moment. This way, the low energy dipole transitions are 
no longer forbidden, and the spontaneous emission rates are gen- 
erally increased. 

In order to investigate the influence of stellar UV radiation on 
the emissivity of Lyman alpha photons, we have used state-of- 
the-art high resolution cosmological simulations and include for 
the first time the detailed chemistry of all relevant processes for 
primordial gas (i.e., photo-ionization, photo-dissociation, colli- 
sional ionization, collisional and radiative recombination and the 
formation of molecules), treating atomic hydrogen as a multi- 
level atom and including the effects of line trapping. In this pa- 
per, we present a number of simulations with different strengths 
of background radiation fields and escape fractions of ionizing 
photons to explore the relative importance of ionizing and photo- 
dissociating radiation. We follow the collapse of a zero metal- 
licity halo and compute the Lyman a emissivity under different 
conditions. 

This paper is structured as follows. In section 2, we briefly 
discuss the simulation setup and numerical methods used. In sec- 
tion 3, we give an overview of our chemical model. In section 4, 
we present the results obtained. Finally, in section 5 we discuss 
our conclusions. 

2. Numerical Methodology 

To investigate the role of background radiation fields in the emis- 
sion of Lyman alpha photons from protogalactic halos, we use 



the code FLASH dDubev et al] l2009h . FLASH is an adaptive 
mesh, parallel simulations code. It is a module-based Eulerian 
grid code which can solve a broad range of astrophysics prob- 
lems. It makes use of the PARAMESH library to handle block 
structured adaptive grids. It uses the message passing interface 
(MPI) library to achieve portability and scalability on many dif- 
ferent parallel systems. The adaptive mesh technique is used to 
achieve high dynamic resolution with minimum computational 
time. We use adaptive mesh refinement (AMR) to add resolution 
in the regions of interest. We use an unsplit hydro solver with 
the 3rd order piece - wise parabolic method (PP M) for hydrody- 
namic calculations (Colella & Woodward 1984). This method is 
well suited for flows involving shocks and contact disc ontinu- 
ities. A multigrid solver is used to solve for self-gravity dRickerl 
120081) . The dark matter is simulated based on the particle mesh 
(PM) method. 

We run the COSMICS package developed by iBertschingert 
(1 19951) to produce Gaussian random field initial conditions. 
We make use of constraint realizations implemented in the 
grafic code (part of the COSMICS package) to select a massive 
halo. We start our simulation at redshift 100 as determined by 
COSMICS. Our computational box has a comoving size of 10 
Mpc. We use periodic boundary conditions both for hydrody- 
namics and gravity. We start the simulations with an initial reso- 
lution of 512 3 grid cells in the central 1 Mpc region and set the 
rest of the box to a resolution of 128 3 grid cells. We initialize 
2.6 x 10 6 particles to compute the evolution of dark matter. For 
the cosmology, we adopt the ACDM model with the WMAP 
5-year parameters (Q m = 0.2581, Ho = 72 km s" 1 Mpc -1 , 
£2b = 0.0441) and a scale-invariant power spectrum with the 
preferred value of erg = 0.8. We select the most massive halo 
and follow its collapse using the AMR method. We impose 7 
additional levels of refinement, which gives 15 levels of refine- 
ment in total. This corresponds to a dynamical resolution limit 
of 75 pc in comoving units. We resolve the Jeans length with a 
minimum of 20 cells. This should be enough to properly resolve 
the essential physics, for further detai l s abou t the significance of 
resolution criteria see lFederrath et al.l (1201 ll) . This automatically 
ensures that we also f ulfill the Truelqve cri terion to suppress ar- 
tificial fragmentation (iTruelove et al.ll 19971) . 

3. Chemistry and Cooling 

Primordial gas is mainly composed of hydrogen and helium. 
Microphysical processes of these atoms play a vital role in the 
formation of the first structures. Stars and galaxies are formed 
in protogalactic halos where collapse is induced due to cooling 
by hydrogen lines. Lyman alpha line cooling is the dominant 
process in metal free gas for temperatures down to 8000 K. H2 
an d HP are the only relevant coolants at temp er atures below 10 4 
K (IGalli & Pallal fl998t iLepp & Shulll [l98l . ISaslaw & Zipovl 
dl967b were first to realize the importance of molecular hydro- 
gen formation through gas phase reactions. They found that trace 
amounts of H2 can be formed through gas-phase reactions al- 
ready at z ~ 300. The rotational and vibrational modes of molec- 
ular hydrogen can be excited at low temperatures and cool the 
gas to a few hundred Kelvin. HD is an important coolant around 
100 K and can cool the gas to lower temperatures. 

We expect that the metallicity in the IGM may still be close 
to primordial at high redshift, in particular for isolated halos. 
For Lyman a emission driven by accretion and infall, adopting 
a primordial composition thus allows us to capture the dom- 
inant effects regulating the emission process. Even for metal- 
licities up to 1 % solar, cooling procedes mostly via primordial 
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Nr. 


UV flux 


Escape fraction 


Lyman alpha Luminosity 


Lyman alpha flux 




J2 1 [erg / cm 2 / s /Hz/ sr] 




[erg Is] 


[erg/cm 2 /s] 


1 


0.01 


0.01 


10 4J 


5 x 10~ 18 


z 


U. 1 


U.Ul 


1 n44 
1U 




3 


10 


0.01 


10 46 


5 x 10" 15 


4 


100 


0.01 


5 x 10 45 


10- 15 


5 


1000 


0.01 


10 45 


2 x 10" 16 


6 


0.01 


0.1 


10 44 


4 x 10" 17 


7 


0.1 


0.1 


10 46 


5 x 10" 15 


8 


0.01 


0.5 


10 46 


5 x 10" 15 


9 


0.1 


0.5 


10 46 


5 x 10" 15 



coolants for the den sities in our simulations (lOmukai et al .12 005: 
IJappsen et al.ll2007l) . Depending on the proper ties of dust grains, 
they c ould potentially enhance H2 formation (ICazaux & SpaansI 
2009). However, we assume here that either the grain-size distri- 
bution is such that this effect is not significant, or that the addi- 
tional H2 formation is suppressed by the photodissociating back- 
ground. The effects of dust will be explored in more detail in an 
upcoming paper. 

To study the thermal evolution of the gas, it is crucial to 
model the chemistry of primordial gas in detail. Solving the 
rate equations self-consistently along with hydrodynamics is a 
computationally demanding task due to the stiff nature of the 
reaction network. We devised a chemical network consisting 
of 36 reactions which involve 23 collisional and 13 radiative 
processes. The reactions and rate coefficients for these pro- 
cesses and their references are listed in table 1 of appendix 
A. We solve the rate equations of the following 12 species 

r, H", H2, h:, d, d + , hd. 



dl997h: 



Abel et al.l 



H, H + , He, He + , He ++ , 1 
Our n e twork is based on lAnninos et al.l . 
(ll997l) : ISchieicher et al.l (l2008bl) : lGlover & Jappsenl (120071) . see 
Appendix A, and includes collisional ionization, radiative re- 
combination, photo-ionization and photo-dissociation rates of all 
the species. The rate equations for these species are solved us- 
ing the backward differencing scheme (BDF), see lAnninos et al.l 
The BDF method is not fully implicit but still more 
accurate and stable than other schemes. We have implemented 
a comprehensive model for cooling and heating of primordial 
gas. It comprises collisional excitation cooling, collisional ion- 
ization cooling, recombination cooling, Bremsstrahlung cooling, 
Compton cooling/heating, photo-ionization heating and photo- 
dissociation heating, as well as a multi-level model for atomic 
hydrogen line cooling. We also included the cooling due to H2 
and HD molecules. The heating / cooling rates are listed in Table 
2 of appendix A. Modeling details of the chemical network are 
summarized in appendix A. 

Another important ingredient of our modeling is the effect 
of line trapping of Lyman alpha photons. Lyman alpha emission 
is an efficient cooling process for halos with virial temperatures 
> 10 4 K. It is oftenly presumed that gas remains effectively opti- 
cally thin to Lyman alpha cooling, but large hydrogen column 
densities produce large optical depths to these lines. At opti- 
cal depths To > 10 7 , the photon escape time typically becomes 
longer than gas free fall time. Because of the weak dependence 
of photon escape time on the gas number density, T p h oc n" 1/9 , as 
compared to free fall time, Tff = 1.0/(Gp ) 1/2 , trapping become s 
effective during the process of collapse dSpaans & Silkl 120061) . 
A detailed description of our treatment for Lyman alpha trap- 



ping and the modeling of the hydrogen line transitions is given 
in Appendix A. 

ISchleicher etaD d2010b) found that in the presence of Ly a 
line trapping, total cooling is compensated as it proceeds through 
the other transitions of atomic hydrogen. The transition from 2p- 
ls produces Lyman alpha photons while the transition from 2s- 
ls produces two continuum photons. These continuum photons 
are not trapped and cooling proceeds through them, depending 
on the ambient density. To model all these mechanisms, we have 
treated atomic hydrogen as multi-level atom by including up to 
the fifth excited state. We have also considered the electronic 
transitions from higher levels of hydrogen. 

The radiation produced during the reionization process can 
have important implications for subsequent structure formation. 
These photons cannot only photo-heat and ionize the gas but also 
can photo-dissociate the H2 and HD molecules. The input spec- 
trum is a blackbody with T* = 10 4 K, with the i ntensity above 
13.6 eV reduced by the escape fraction f esc as (lOmukai et al] 
120081) 



= J 2 il0- 21 [B y (T,)/B yH ]f(v) erg cnrW'Hz 



-2-1 -li 



.-1 



(4) 



where J21 is the mean intensity of radiation at the Lyman 
limit, T* is the color temperature of a star. f(v) is one for non- 
ionizing photons and equal to f esc for ionizing photons. B y (T) is 
the stellar radiation intensity (i.e. black body spectrum), 



B V (T) 



2hv 3 



1 



c 2 (exp(hv/kT)- 1) 



(5) 



where h is Planck constant, v is the frequency of radiation, 
k is Boltzmann constant and T is the radiation temperature. 
Jih = hi x fesc represents the ionizing soft UV flux and J21 is 
the background UV flux below the Lyman limit (i.e. 13.6 eV) 
in units of 10" 21 erg /cm 2 /s/Hz/ sr. While propagating through 
the intergalactic medium, the photons are attenuated as the op- 
tical depth depends on the frequency. In order to approximately 
account for the effect of self- shielding of gas in the halo, we es- 
timate the optical depth as a function of frequency as 



Ti(vj) = licri(vj) X ni X Aj 



(6) 



where i denotes the species, j the frequency bin, cr the photo- 
ionization cross- section of the given species, ni its number den- 
sity and A] is the local Jeans length. For a given optical depth, 
attenuation of the radiation field is calculated as 



J att (y) = r xt (y)xexp(-r(y 1 )). 



(7) 
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Fig. 1. The panels in this figure show the density and temperature slices at redshift 5.5. The left-hand panels show density slices 
through the center of the halo for J21 = 1000, 10, 0.1 from top to bottom and the right-hand panels show the corresponding 
temperature slices. Velocity vectors are overplotted on the density slices. The figure shows the inner region of 20 kpc in comoving 
units. All plots shown here are for an escape fraction of 1%. 
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Fig. 2. Figure shows the density radial profiles for different escape fractions. The left panel shows density profiles for an escape 
fraction of 1%. The right panel shows the density radial profile for escape fractions of 10% and 50%. 
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Fig. 3. Left panel shows the radial velocity profile for different background UV radiation fields as given in the legend. Right panel 
shows the tangential velocity radial profile of the halo for different background UV fields. For the tangential velocity, the average 
was calculated over the absolute values. So, it is a measure of turbulence rather than rotation. 



4. Results 

We perform nine cosmological simulations with different 
strengths of the background UV field and escape fractions of 
ionization radiation as listed in table [T] We start our simulations 
at redshift 100. We note that density perturbations decouple from 
the Hubble flow and begin to collapse through gravitational in- 
stabilities. During the non-linear evolution phase gas collapses 
in the dark matter potentials and gets shock heated. We select 
a halo of 10 9 M© at redshift 10 and follow its collapse down 
to redshift 5.5. At a redshift of about 8, gas becomes photo- 
ionized as we switch on the ionizing background produced by 
reionization. The temperature of the gas is raised above 10 4 K 
due to photo-ionization heating depending upon the strength of 
the background flux. Photo-ionized gas recombines and begins 
to cool and collapse because of recombination cooling. The viri- 
alization process continuously transforms the gravitational po- 
tential energy into kinetic energy of gas and dark matter. Due to 
the dissipative nature of the gas, part of the gravitational energy 
goes into thermal energy and is radiated away. During virializa- 
tion, the gas becomes turbulent and settles in the center of a halo 
through cold streams as shown in figure [T] The temperature of 
the gas in filamentary accretion is ~ 10 4 K and densities are 
10" 2 - 1 cm -3 . The density radial profiles are shown in figure [2] 
It can be seen from the figure that density varies with distance 
r like r -2 8 . Our simulations reach a maximum density of a few 



times 10 5 cm -3 . Small fluctuations in the density radial profiles 
are due to the dumpiness of the medium as the halo is not per- 
fectly symmetric. We show the radial velocity profile for three 
selected cases in the left panel of figure [3] The figure shows that 
gas is falling into the center of the halo. The variations in the 
radial velocity profile at the outer radii are due to different ac- 
cretion rates for various background radiation fields. The gas in 
the center of the halo becomes stable. The right panel of figure [3] 
shows the tangential velocity radial profile. Variations in tangen- 
tial velocity are caused by shocks. The total rms velocity hardly 
depends on the radiation field, as it is mostly due to the virial 
temperature. We find typical velocities of about -50 km/s as 
expected from the virial temperature. Virialization shocks also 
play an import ant role in the emis sion of Lyman alpha, for fur- 
ther details see lLatif et all (l2011ah . 

We see that for J21 < 10, the temperature of the gas is ~ 10 4 
K while for higher intensities (i.e. J 21 > 10) it reaches 10 5 K. 
The temperature of 10 5 K is driven by the ionization of helium, 
which absorbs more energetic photons. Cooling due to helium 
lines becomes effective at temperatures above 5x 10 4 K and keeps 
the gas temperature at ~ 10 5 K. At densities of 1 cm -3 , cooling 
due to radiative recombination becomes very effective and cools 
the gas down to 10 4 K. The gas in the surroundings of the halo 
remains hot (> 10 4 K) as shown in the top left panel of figure 
ID Similarly, the degree of ionization varies depending upon the 
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Fig. 4. The escape fraction of ionizing radiation is 1% for all the panels shown. Colors show different values of background 
intensities J21 as indicated in the legend. The upper left panel of this figure shows the temperature radial profile of the halo. The HII 
abundance radial profile for the halo is depicted in the upper right panel. H2 abundance is shown in the lower left panel. The Hell 
radial profile of the halo is depicted in the lower right panel. It can be noted that for higher values of J21, the IGM gas is heated up 
to 10 5 K and the degree of ionization is enhanced. H2 is photo-dissociated for stronger radiation fields. 



impinging radiation intensity, as shown in the upper right panel 
of figured It is small for lower fluxes and large for higher fluxes. 
We also see that helium remains mostly neutral for lower fluxes 
and becomes singly ionized for a stronger radiation background. 
The He + abundance is shown in the bottom right panel of figure 
[4] At densities > 1 cm -3 , the recombination rate becomes high 
and the gas begins to recombine. Consequently, the degree of 
ionization declines. 



For low background fluxes, a H2 abundance of 10" 6 is 
formed at low densities (i.e. 10" 2 cm -3 ) and a temperature of 10 4 
K, while its formation remains inhibited for higher fluxes at low 
densities. For a high background radiation field, i.e., J21 > 1000, 
H2 is photo-dissociated even in the center of the halo. For weaker 
radiation fields photo-dissociation of molecular hydrogen be- 
comes ineffective and the H2 abundance reaches the universal 
value (i.e., 10" 3 ). The H2 radial profile for different fluxes is 
shown in the bottom left pane l of figure |4] Our results are con- 
sistent with previous studies (iShang et al. 1201 (jt iDijkstra et al.l 
120081: iBromm & Loebll200l I Johnson et al.ll201ll) . We also com- 
puted the column densities for different background fluxes, find- 
ing that the column density is higher for low fluxes in the enve- 
lope of a halo due to a higher neutral gas fraction. Conversely, 
the column is low for higher radiation field strengths because of 
a higher ionization degree. The upper left panel of figure[5]shows 
the column density radial profile. 



Based on these results, we computed the emissivity of 
Lyman alpha photons, finding that it is enhanced in the pres- 
ence of a background incoming radiation field. It increases with 
the strength of the radiation field up to 10" 22 erg/cm 2 /s/Hz/sr 
(i.e. J21 x f esc ) and then starts to decline for higher intensity. This 
decline in emissivity is a consequence of the decreasing abun- 
dance of atomic hydrogen. At columns higher than 10 22 cm -2 , 
trapping of Lyman alpha becomes effective and consequently the 
emission of Lyman alpha photons diminishes while cooling still 
proceeds through higher electronic transitions of hydrogen. For a 
fixed ionizing background, the emission increases with increas- 
ing photo-dissociation flux. The emissivity of Lyman alpha pho- 
tons as well as the total emissivity is plotted against radius in the 
upper right panel of Fig. [5] Lyman alpha emission is higher in 
the envelope of the halo and sharply declines towards the center. 
This behavior of Lyman alpha emission is due to line trapping 
of Lyman alpha photons and is consistent with previous stud- 
ies (lLatif et al.l l2011at [Schleicher et all l2010bt ISpaans & Silkl 
l2006ir In addition to the line trapping, also the cooling by H2 
and HD may prevent Lyman alpha emission, as it drives the gas 
to lower temperatures (i.e. the run J21 = 10). The total emissiv- 
ity remains roughly constant as cooling is compensated by two- 
photon continuum processes. Moreover, for higher background 
radiation field strengths the total emissivity from hydrogen elec- 
tronic states is higher than for low radiation fields. This is due to 
a higher molecular hydrogen fraction for low fluxes and a lower 
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Fig. 5. The upper left panel shows the column density radial profile of the halo. The emissivity radial profile for the halo is depicted 
in the upper right panel, where the dot-dashed line shows the Lyman alpha emissivity and the triple dot-dashed line shows the total 
emissivity of all the coolants. The lower left panel shows the luminosity radial profile. The flux radial profile of the halo is depicted 
in the lower right panel. All panels in this figure have an escape fraction of 1% for ionizing photons. Colors represent the different 
values of background intensity J21 as shown in the legend. Emissivity of Lyman alpha photons increases upto J 21 = 10 and starts to 
decrease for intenser radiation fields as the neutral fraction goes down. Lyman alpha luminosity and flux show the same behavior. 



H2 fraction for high fluxes through photodissociation. The lo- 
cal variations in emissivity are because of the dumpiness of the 
medium. 

The enclosed Lyman alpha luminosity as a function of ra- 
dius is shown in the bottom left panel of figure \5\ The lumi- 
nosity increases with the strength of the ionizing radiation field 
up to 0.1 (in units of J21) and then decreases with increasing 
radiation field strength similar to the emissivity. Luminosity val- 
ues are of the order of 10 43 - 10 46 erg/s which are compara- 
ble to observed Lyman alpha blobs as w e ll as estimates from 
numerical simulat i ons (lOuchi et al.l 120091: iLehnert et al.1 1201 Ot 
iGoerdt et afll201(t IPiTkstra et al.ll2006ah . We also computed the 
flux. We found that the flux mostly emerges from the envelope of 
a halo rather than its center and is significantly enhanced in the 
presence of a background UV field. The total flux varies from 
10" 18 - 10" 15 erg/cm 2 /s depending upon the ambient radiation 
field strength. Our results are in agreement with our previous 
finding that emission of Lyman alpha is extended and emerges 
from the hot gas in the surroundings of a halo. 

Due t o large uncertainties in th e escape fraction of ionizing 
radiation (iFernandez & ShulfeOlOl) . we have performed a num- 
ber of simulations to assess the impact of different escape frac- 
tions on the emissivity of Lyman alpha photons. We performed 
simulations with escape fractions of 1%, 10% and 50 % for a 
background radiation field of J21 = 0.01 and 0.1. We found 



that for a similar value of J21 describing the strength of pho- 
todissociation, an increasing escape fraction of ionizing photons 
enhances the temperature of the gas. Radial profiles of temper- 
ature, H + , He + and H2 abundances are shown in figure [6] As 
expected, the degree of ionization increases for higher escape 
fraction. The gas is still neutral in the center of the halo but be- 
comes ionized in the surroundings. Column density is shown in 
the upper right panel of figure [71 It is similar to the case with 
low escape fraction. Differences in the columns is due to dif- 
ferent gas ionization fractions. The upper left panel of figure [7] 
shows the emissivity of Lyman alpha photons. It is found that for 
the same background value of J21 , the emissivity is enhanced for 
higher escape fractions. It reaches a maximum in the envelope 
of the halo and sharply drops as trapping becomes effective in 
the center. The total emissivity is again compensated by cooling 
through higher states of hydrogen atoms. There are some local 
variations in the emissivity profile due to the density structure, 
fluctuations in temperature and the ionization degree. We also 
computed the luminosity and flux for this case and found that 
the flux is boosted for a higher escape fraction. The radial pro- 
files of luminosity and flux are shown in the bottom panels of 
figure [7] Luminosities and fluxes for all cases are summarized in 
tableffl 

We also computed the emissivity of various coolants to 
assess their relative importance for the cooling. We here 
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Fig. 6. The upper left panel of this figure shows the temperature radial profile of the halo. The HII abundance radial profile for the 
halo is depicted in the upper right panel. H2 abundance is shown in the lower left panel. The Hell radial profile of the halo is depicted 
in the lower right panel. The legend in each panel shows the values of the escape fraction of ionizing photons. The corresponding 
strength of the background radiation field is also shown in the legend. It is found that using a higher escape fraction of ionizing 
radiation, the degree of ionization is enhanced and gas is heated upto higher temperatures for the same background radiation field 
strengths. The formation of H2 is not affected much by a higher escape fraction of ionizing photons. 



show three typical cases for background UV fluxes of J 21 = 
0.01, 0.1 and 1000. Radial profiles of cooling rates from H(2p - 
Is), H(2s - Is) and H(3 - 2) transitions as well as H2 and HD 
ro- vibrational transitions are shown in figures [8] and 13 The left 
panel of figure [8] shows the cooling rates for J21 = 1000. It 
is found that at radii > 3 x 10 20 cm cooling is only coming 
from electronic transitions of hydrogen. Notable transitions are 
H(2p-ls), H(2s-ls) and H(3-2). The population of these transi- 
tions is regulated by radiative de-excitation and collisional ex- 
citation, which makes them efficient coolants. The cooling from 
H2 and HD remains negligible at these radii. At radii between 
3 x 10 20 - 10 20 cm line trapping of Lyman alpha photons becomes 
effective and two-photon continuum emission (2s- Is transition) 
becomes dominant. In the center of the halo (< 10 20 cm) cool- 
ing is regulated by trace amounts of molecular hydrogen. The 
right panel of figure [8] shows the cooling rates for J21 = 10. At 
radii > 10 21 cm cooling is dominated by electronic transitions of 
hydrogen. Molecular hydrogen becomes an effective coolant for 
radii between 10 21 - 10 20 cm as cooling from electronic transi- 
tions of hydrogen becomes negligible due to lower gas temper- 
atures. In the core of the halo, HD cooling dominates over H2 
cooling for J21 < 0.1. It can be seen in figure [9] that, similar to 
previous cases, at radii > 10 21 cm the halo is cooled by elec- 
tronic transitions of hydrogen. Emission from these transitions 
becomes negligible at lower radii due to lower gas temperatures 



(< 10 3 K). Again, cooling in the center of the halo is dominated 
byHD. 

5. Discussion and Conclusions 

We have performed 9 simulations in total to investigate the 
role of the UV background radiation field in the emission of 
Lyman alpha photons. We used the adaptive mesh refinement 
code FLASH to carry out cosmological simulations and cou- 
pled it with a chemical network that solves the rate equations 
of 12 species. Our chemical network includes collisional ion- 
ization, radiative recombination, photo-ionization and photo- 
dissociation processes, the formation of molecules and a multi- 
level model of atomic hydrogen. We found that the presence of 
a background flux photo-ionizes the gas and raises the gas tem- 
perature to > 10 4 K. We found that for weak radiation fields 
(J21 < 100), cooling due to H2 decreases the temperature to at 
least 1000 K. For J21 < 0.1, cooling by HD is important as well, 
leading to gas temperatures of ~ 100 K . We note that t he HP 
cooling function was recently revisited by Coppol a et al ] d20TTh . 
However, they only find differences at high temperatures, where 
H2 dominates the cooling. The thermal evolution in these cases is 
remarkably different and provides an environment conducive to 
star formation in the halos irradiated by weaker radiation fields, 
while stronger fields may lead to the formation of massive ob- 
jects. 
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Fig. 7. The upper left panel shows the column density radial profile of the halo. The emissivity radial profile for the halo is depicted 
in the upper right panel, where the dot-dashed line shows the Lyman alpha emissivity while the triple dot-dashed line shows the 
total emissivity of all the coolants. The lower left panel shows the luminosity radial profile. The radial profile of the enclosed flux 
from the halo is depicted in the lower right panel. The escape fraction of ionizing radiation and the strength of background radiation 
field is shown in the legend of each panel. The emissivity of Lyman alpha is significantly enhanced for higher escape fraction and a 
given background UV field strength. The same holds for luminosities and fluxes. 
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Fig. 8. Cooling rates of H2, HD, as well as H(2p-ls), H(2s-ls) and H(3-2) transitions are plotted against radius in both panels. Left 
panel is for J21 = 1000 with a 1% escape fraction of ionizing radiation. Right panel is for J21 =0.1 with an escape fraction of 10%. 



We noted that the ubiquity of background flux enhances 
the emission of Lyman alpha photons. We found that the emis- 
sion of Lyman alpha photons is boosted with ionizing flux up 
to 10" 22 erg/cm 2 / s/Hz/ sr. The emission starts to decrease for 
higher values of the ionizing flux. This decline in the emissivity 
of Lyman alpha photons is because of the small abundance of 
atomic hydrogen in the presence of intense radiation fields. The 



flux predominantly arises from collisional excitation of atomic 
hydrogen. We found that the contribution from the recombina- 
tion is at least few orders of magnitudes smaller. 

We also noted that the emission of Lyman alpha photons 
emerges from the envelope of a halo rather than its center. At 
columns above 10 22 cm -2 trapping of Lyman alpha photons be- 
comes effective and consequently the emissivity drops sharply 
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above t hese columns. Our results are consistent wi t h previous 
studies (ISpaans & Sil3l2006l:[Schieicher et al.ll2010bl:lLatif eTaD 
l2011bl) . We also computed the flux emanating from the halo. We 
find that the flux is of the order of 10" 17 - 5 x 10" 15 erg/cm 2 /s 
for different background UV fluxes. Similar to the emissivity, 
the flux is extended and comes from the relatively hot gas in 
the outskirts of a halo. The Lyman alpha wavelength observed 
at redshift 5.5 will be 0.79 /im. This should be detectable with 
the upcoming telescope JWST. NIRcam/NIRSpec instruments 
available on JWST will be well suited for its detections. The 
mass of our halo is 10 10 M© and has an angular size of 0.52 arc- 
sec at redshift 5.5, attainable by JWST. Subaru telescope will 
also be able to detect this flux using suprime-cam in combination 
with intermediate-ba nd filters (IA679, IA 505) as discussed in the 
previous studies, see ISaito et all (120061) . Our results are in the 
range of planned surveys in future with Subaru which will make 
use of next-generation wide-field camera, Hyper Suprime-Cam, 
and several narrowband filters. We also computed the surface 
brightness profiles for different background UV fluxes which are 
shown in figure \TU[ The surface brightness profiles are flat as 
most of the flux emerges from the envelope of a halo rather than 
its core. They are comparable to the previously observed values 
(ISaito et al.ll2006l: lOuchi et al.ll2Q09h . 

In the presence of large velocity gradients, line trapping of 
Lyman alpha photons does not remain effective. For r > 10 7 , 
only photons far-away in the line wings will be able to escape. 
Trapping in our case occurs at radii > 10 205 cm, where large ve- 
locity gradients are significantly reduced. Scattering of Lyman 
alpha photons which perform a random walk in frequency and 
space space will blue shift the frequencies l eading to an up- 
per li mit on the trapping time ~ tff(Av/v) dRees & Ostrikerl 
fl977h . Although, there are some uncertainties in the diffusion 
time scale. The issue of trapping has been studied in detail by 
ISpaans & Silkl (120061) . They found that weak dependence of pho- 
ton escape time on number density (t esc n" 1/9 ) as compared 
to the free fall time results in t esc /tff ~ n 7/18 . Consequently, trap- 
ping will remain effective. 

We found that the critical value of the background UV 
flux for H2 formation is J21 = 1000, which is in accordance 
with previous estimates dOmukaill200"ll; iBromm & Loeb 2003; 
IWise & Abelll2007l:IShang et alJl2O10h 7 We also noticed that the 
HD abundance remains lower than the critical value (i.e., 10" 6 ) 
for J21 > 0. 1 . It could be due to little amount of H2 formation and 
direct photodissociation of HD. Consequently, HD cooling be- 



comes important for lower values of background flux (i.e., J21 < 
0.1). Our results ar e consistent with the critical value of UV flux 
for HD derived by lWolcott-Green & Haimanl (1201 ll) . Collisions 
produce trace amounts of HD at high densities which is too 
small to dominate the cooling for higher values of J21 . It has 
been found that HD cooling is mainly controlled by H2 forma- 
tion rather than direct dissociation. For further details on the role 
of UV radiation in the formation of HD and its implications for 
the str u cture formation see dedicat e d stud ies iMcGreer & Bryanl 
(120081) : IWolcott-Green & Haimanl (120111) . We also performed 
simulations to take into account the local variations in the back- 
ground flux, i.e., J21 = 100 and 1000. We found that for these 
cases, emission of Lyman alpha photons is reduced, even for an 
escape fraction of even 1%, compared to cases with lower val- 
ues of J21 =0.1 and 0.01. However, these are rare occurences, 
which will only apply in the near vicinity of sources (immedi- 
ate neighbours). We also studied a few cases for different escape 
fractions up to 50%, and found that a higher escape fraction of 
ionizing photons enhances the emissivity of Lyman alpha pho- 
tons. We also computed the emission in the helium Lyman alpha 
line and found that it is considerably weaker than HI Ly a, due 
to the low abundance of ionized helium in the halo. This may, 
however, be different in the presence of an AGN or HII regions 
from m assive stars. Furthe r di scussion on helium l ines can be 
found in I Yang et al] (120061) and lHaiman etall d2000h . 

In our case, ionization predominantly occurs at low densi- 
ties. Therefore, the contribution of the Lyman alpha recombina- 
tion emission is a few orders of magnitude lower than the col- 
lisionally excited line emission. At radii 10 205 cm, emission di- 
minishes as Lyman alpha line trapping becomes effective in the 
center. In the case of radiation sources at high densities, recom- 
bination emission could be more important. We have further as- 
sumed that the background UV flux is constant in time, although 
we do not expect significant differences if a more detailed model 
of the time evolution of the UV flux is adopted. iDijkstra et al.l 
(12004]) found that the impact of a radiation background is less 
at high redshift (z > 10) as compared to z = 3 due to different 
properties of the halo. 

In light of these results, it will be important to obtain a better 
understanding of the environmental conditions, in particular the 
strength of the UV background radiation field, for ~ 10 10 M© 
halos close to the reionization epoch. We are further planning to 
assess additional uncertainties, in particular those related to the 
metallicity of the gas. 

We have assumed here that the halo is metal free for 
the reasons mentioned in sect ion 3. Indeed pr istine halos 
may exist down to redshift 6 (iTrenti et al.l 120091) . The pres- 
ence of small amounts of dust may absorb the Lyman al- 
pha photons and re-emit them as continuum at lower fre- 
quencies. This can change the scenario. Nevertheless, if the 
medium is i nhomogeneous, escape of Lyman a bec omes pos- 
sible again { Haiman & Spaans 1999; Neufeld 19 9l|). also see 
iGould & Weinberg! (ll996h : lHansen & Ol3 (l2006h : IShimizu et all 
(1201 ll) . H2 formation is significan tly enhanced even for low dust 
abundances of 10" 5 - 10" 4 solar dCazaux & Spaansll2009l) . The 
impact of dust on Lyman alpha emission is an important issue 
and needs a dedicated study. In an upcoming paper we will ex- 
plore the influence of dust on the emission of Lyman alpha in de- 
tail. The ubiquity of outflows may even enhance the emission of 
Lyman alpha ph otons up to 5%, by directl y injecting them along 
the line of sight (iDijkstra & Wyithell2010l) while the presence of 
turbulence may amplify magnetic fields in the se protogalaxies, 
which may affect subsequent star formation (iSchleicher et al.l 
I2010al) . . 
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Fig. 10. Both panels of this figure show the surface brightness profiles for various strengths of background UV radiation fields as 
mentioned in the legends. The value of the escape fraction is 1% for the left panel. For the right panel, the value of the escape 
fraction against each radiation field is shown in the legend. 



We adopted a background UV flux with a T* = 10 4 K ther- 
mal spectrum, typical for radiation dominated by Pop II stars 
when z < 6. The H2 abundance is mainly controlled by H" pho- 
todissociation for a soft UV spectrum in our simulations. For 
higher radiation temperatures, it will be governed by H2 pho- 
todissociation. This may have a strong impact on the critical J21 
values of H2 and HD. If Pop III stars with T* = 10 5 K domi- 
nate the UV background, then the critical value of J21 required 
for the photo-dissociation of H2 increase s by two orders of mag- 
nitude as computed by previous studi es (iBromm & Lo eb 2003; 
IWise & Abelll2QQ7fc IShang et al.ll20Tol) . On the other hand, for a 
T* = 10 5 K blackbody spectrum, an almost 100 times smaller 
value for the ionizing UV background suffices to photo-ionize 
the gas. So, our estimates may vary for a harder spectrum (i.e., 
10 5 K) but we do not expect a big change in the overall conclu- 
sions. 

Our results include local radiative transfer effects but 
nonetheless are comparable to c omplete radiative trans- 
fer studies (Diikstra et al. 1 l2006allbl: iThomas & Zaroubil [2 008: 
Wise & Abel 2008b; Verha mme et alJ 120081: iMaselli et all 
200 91 iHocuk & Spaansl [20101: iFaucher-Giguere et al.l l2010b 
Zheng etal.1 120101 l201ll) . Lyman alpha radiative transfer 
strongly depends on the density and velocity profile of a halo due 
to the resonance nature of the Lyman alpha line and its large scat- 
tering cross section. Therefore, we expect some variation in the 
emergent Lyman alpha flux due to the presence of anisotropies 
in HI columns along different lines of sight. 

We have also ignored the impact of X-rays on the emission of 
Lyman alpha photons. The presence of X-rays will heat the gas 
and increase the emissivity of Lyman alpha. Lyman alpha blobs 
powered by X-ray sources ha ve already been det ected (roughly 
17% of the observed sources, iGeach et aD l2009h . Due to com- 
plex chemistry involved in the modeling, X-ray sources demand 
complete radiative transfer simulations to probe their impact on 
Lyman alpha emission. We intend to pursue this in the future. 
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Appendix A: 

A.1. Chemical Model 

We self-consistently solve the chemical rate equations 
for our cosmological simulations in an expanding uni- 
verse. Rate equations are solved for the following species: 



H, H\ He, He + , He + 



H" 



H 2 , H+, D, D + , HD under 



non-equilibrium ionization/equilibrium conditions. We also 
treat electronic states of hydrogen, up to 5 levels, as separate 
species. The rate equations are 



diii _ _ 
— - = -Dmi + Q, 
dt 



(A.l) 



where ni is the number density of the ith specie, and Di and Q are 
destruction and creation coefficients for the ith specie. Equation 
lA.ll is solved using the backward differencing scheme (BDF). 
The solution for equation lA.ll is 



C new dt + n: 



old 



ni = 



1 + D new dt 



(A.2) 



Values for C^ ew and Dj* ew are approximated using species den- 
sities from the previous time step and those from the new time 
step. During each hydrodynamical time step, the chemical time 
step is computed as 



dt c = min(HI/(dHI/dt),H + /(dH + /dt))/10. 



(A.3) 



For a given chemical time step, species are evolved until 
the chemical time step becomes equal to the hydrodynami- 
cal time step. We take the initial abundances of species from 
ISchleicher et al] (l2008bl) . These are consistent with nucleosyn- 
thesis and observational constraints. 

Strong radiation fluxes can photo-ionize the gas, and for 
a primordial composition specifically hydrogen and helium. 
Photo-io nization rates are cal culated by solving the following 
equation ( A nninos et al . 199^): 



dp 
di 



Ipk 



Jvth 



hv 



(A.4) 



where k denotes H, He, He + , I(v) is the intensity of the radi- 
ation field, v t h is the threshold energy for each rate for which 
ionization takes place and cr\ is the frequency dependent photo- 
ionization rate. In order to compute the photo-heating produced 
due to ionization of the gas the following equation is solved: 



T= f 4;rcr^(hv-hv th )dv. 



(A.5) 



Details about H2 and HD chemistry are given in the section be- 
low. 



A.2. H 2 and HD chemistry 

Neutral hydrogen can capture free electrons, which leads to H2 
through the H" route. Another channel for H2 formation is the 
collision of H atoms with H+. Similarly, HD can be formed 
through collisions of H2 molecules with deuterium species. We 
have included collisional and radiative formation and destruc- 
tion of H2 and HD in our network. The relevant rates are listed 
in table 1 . We have also incorporated the photodissociation rates 
of H", H2 and H+ by the background radiation field. Rate equa- 
tions for H2 and HD species are solved under non-equilibrium 
conditions, while for the ions H" and H+ chemical equilibrium is 
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established rapidly. The reactions involving these species occur 
on short time scale. The photodissociation and photo-ionization 
rates of these species are calculated by solving equation IA.41 
where k denotes H, H", and H2, v t h is the threshold energy 
for photo-dissociation and a\ is the frequency dependent pho- 
todissociation rate. We also include heating due to photodisso- 
ciation (for both the Solomon process and direct dissociation) 
of H2 by absorption of back ground radiation in the Lyman- 
Werner band. Similarly, we compute the photo-dissociation of 
HD molecules in the presence of a UV field. We take into ac- 
count the self- shielding of H2 in the Lyman- Werner band when 
computing the H2 phot odissociation rate. We use the expression 
given m equation lA31(Sh ang et al.ll201Qh . We also include the 
photodissociation of HD, which is approximately 10% higher 
than the H2 photodissociation rate for the Solomon process. For 
direct dissociation they are roughly the same . To compute the H2 
column density we use the expression from lShang et al] (1201 Ol) 
as given in equation \Kj\ 



fsh = min[l,(- 



N H2 



10 14 cm" 2 

Nh 2 = fH 2 n tot^J, 



-3/4-1 



(A.6) 
(A.7) 



the gas, and (5 is a geometrical factor with J3=3 for spheri- 
cal collapse. Note that t P h ^ Wi for most transitions. For 
the Lyman alpha line and other direct transitions to the ground 
state s expressions fo r optic al depth and diffusion time are taken 
from ISpaans & Silkl (120061) . For Lyman a, t P h = L(ar 2 i) 1/3 /c, 
where a is natural-to-thermal line width and the optical depth 
r 2 i = 1.04 x 10- 13 N H T-°- 5 . The cooling function for hy- 
drogen lines is computed from the escape probabilities, level 



Further details of all these 



> can be found in 


Omukail 


; Rees & Ostrikei 


(I1977D: 



(2010b). 



where Nh 2 is the H2 column density, n tot is the total particle den- 
sity and Aj is the Jeans length, is the self- shielding fa ctor to 
correct the impinging UV flux (Draine & Bertoldi 1996). This 
allows us to compute Nh 2 and f S h from local quantities only. For 
m olecular hydrogen co oling we have used the cooling function 
of iGalli & Pallal (U998I) . We also include he ating/cooling rates 
due to formation and destruction of H2 from lGlover & Jappsenl 
( 2007). We also include HD photodissociation by a background 
UV flux in our c hemical model. The HD p hotodissociation rate 
is adopted from iGlover & Jappsenl (|2007[). For HD cooling w e 
have used a fit to the cooling function of llipovka et al] (120051) . 



A. 3. Lya Trapping 

We have modeled the hydrogen atom as a multi-level sys- 
tem with five levels. Ou r level population model is based on 
ISchleicher et al.l (l2010bl) . For the first excited state we distin- 
guish the levels 2s and 2p. Transitions from 2p-ls produce 
Lyman alpha photons, while transitions from 2s- Is produces the 
2-photon continuum. We calculate transition rates from level i to 
level j using the equation: 

Rij = Ay&se.ijO + Qij) + Qji > j, (A.8) 

Rij = ^Aji&scjiQji + Qj, i < j, (A.9) 

where Ay are the radiative decay rates, /3 e sc,ij the escape proba- 
bilities of for transition ij and C j, the coll i sional excitation and 
de-excitation rates, adopted from lOmukail (1200 ll) . gi the statisti- 
cal weight of the ith level and Qji = c 2 J con t,ij/(2hv?), where v y is 
the frequency of the transition and J con t,ij is the average intensity 
of the background radiation field. We take the b ackground radia- 
tion fi eld to be a black body temperature 10 4 K (iGlover & Savinl 
120091) . The escape probability for the Ly a photons is given by 

1 - exp(-Tjj) -£W 
A>sc,ij = exp(- ), (A. 10) 



where T y - is the optical depth at line center of transition ij, 
t P h is the photon diffusion time, t co n is the collapse time of 



14 



Latif et al.: Lyman alpha emission from the first galaxies 



Table A.l. 



Collisional and radiative rates 



No. Reaction 



Rate coefficient cnrs /Cross-section cr(cm 2 ) 



Reference 



9 
10 
11 
12 

13 
14 
15 
16 
17 
18 
19 
20 
21 

22 
23 

24 
25 
26 

27 

28 
29 



30 



31 

32 
33 
34 
35 
36 
37 



H + e"^H + +2e" 
H + +e"^H + y 
He + e"^ He + + 2e" 
He + + e~— > He + y 

He + + e"-> He ++ + 2e" 
He ++ + e"^He + +y 
H + e-^H" +y 



H" + H^ H 2 + e" 



H + H + ^H 



H + HJ- 



H" 



+ H + 
+ e~- 



2 +r 

H 2 +H + 
>H + H 
H + H 



HJ +H 



H 2 +H + 
H 2 + e~— > H + H + e" 
H 2 + H^ H + H + H 
H 2 + H 2 ^ H 2 + H + H 
D + e"-> D + + e" + e" 
D + +e~^D + y 
H~ + e~-> H + e~ + e~ 
H" + H^ H + H + e" 
H" + H— > H 2 + e" 



H + D + - 
D + H + - 



D + H + 
H + D + 



H 2 +D + ^HD + H + 
HD + H + ^H 2 +D + 
H 2 + HD + H 

HD + H^ D + H 2 



H" +y 
H 2 + y " 



► H + e" 

h: + e" 



HJ + y -> H + H + 



HJ + y -> 2xH + e" 

H 2 + y -> H 2 * -> H + H 
H 2 + y -> H + H 

H + y -> H + + e" 
He + + y -> He ++ + e" 
He + y -> He + + e~ 
HD + y -> H + D 



see reference 
see reference 
see reference 

3.92 x 10" 13 /T^ 6353 for T, 
1.54 x 10" 9 T- L5 exp( 
see reference 

-IOt-0.5/ 



< 8000 

44^9eV )[0 .3 



+ exp(M^)]forT e > 8000 



3.36x10-^(^)^(1 +( T ^r 1 
dex[-17.845 + 0.762(log T) + 0.1523(log T) 2 - 0.03274(log T) 3 ] for T < 6000 
dex[-16.42 + 0.1998(logT) 2 - 5.447 x 10" 3 (logT) 4 + 4.0415 
xl0- 5 (logr) 6 ]for T> 6000 
1.5xl0" 9 forT<300 
4.0 x 10- 9 T-° 17 for T < 300 for T > 300 
dex[-19.38 - 1.523(logT) + 1.118(logT) 2 - 0.1269(logT) 3 ] 
6.4 x 10" 10 

5.4 x 10- 6 T-° 5 + 6.3 x 10" 8 - 9.2 x lO^T 05 + 4.4 x 10" 13 T 

1.0xl0" 8 for T<617 

1.32 x 10- 6 T-°- 76 for T > 617 

see reference 



3.73 x 10- 9 T° 



*exp(- 



6.67xlO- 12 T°- 5 exp(-l +^22) 
_exp(-l + ^) 



[88)0.407^-2.242 



(1.0+6.761xl0- 6 T) 5 - 688 

see reference 

2.753 x 10- 14 (^f^) L5 ((1.0- 
see reference 
see reference 

6.9 xlO- 9 T-° 35 for T< 8000 
9.6 xlO- 7 T-°- 90 for T> 8000 
2.06 x 10- 10 T a396 exp(^) + 2.03 x 10 

2.0 x 10- 10 T a402 exp(^) - 3.31 x 10" 
3.44 x 10- 10 T - 35 for T > 2.0 x 10 5 
[0.417 + 0.864 log T - 0.137(logT) 2 ] x 10 

1.1 x 10- 9 exp(^p) 
1.69 x 10- 10 exp(4P) for T < 200 

-10^~/-4680 , 198800 



-9-p-0.332 

17 T 1.448 forT < 2.0X10 5 



1.69 x 10- lu exp(^ + ™) for T > 200 
5.25 x 10- n exp(^) for T < 200 
5.25 x 10- n exp(^P + ^) for T > 200 
o- 2 3 = 7.928 x 10 5 (y- v th )i(^cm 2 for hv > hy th = 0.775 eV 
o- 2 4 = for hv< 15.42 eV 

o- 24 = 6.2 x 10" 18 hv - 9.4 x 10" 17 cm 2 for 15.42 < hv < 16.50 eV 
o- 24 = 1.4 x 10" 18 hv - 1.48 x 10" 17 cm 2 for 16.50 < hv < 17.7 eV 
o- 24 = for hv > 17.7 eVcm 2 for hv > 17.7 eV 
o-25 - dex[-40.97 + 15.9795(^-) - 3.53934(^) 2 + 0.2581155(^) 3 ] 

1 th 1 th 1 th 

for Eth = 2.65 eV2.65 < E < 11.27 eV 

dex[-30.26 + 7.3935(^) - 1.2914(^) 2 + 6.5785 x 10" 2 (^) 3 ] 
for 11.27 <E< 21.0 eV 
AAZN97 

log(o- 26 ) = -16.926 - 4.528 x 10 2 hv + 2.238 x 10" 4 (hv) 2 
for 30 eV < hv < 90 eV 

1.12 x 10 8 J(v) where J the is flux at E= 12.87 eV 



■ 4.245 x 10- 7 (hv) 3 



C~28 
CT20 



+ o-^°)1.0 — ^- (o-^g 1 + tf - ^ 1 ) for sigma values see reference 



J-(rr U) 

y+ lV^28 1 ^28 
Ao / _v_ \4 exp[4.0-4.0(arctan6)/ , 6] 
Z 2 ^v t j/ l-exp(-2;r/e) 
— A / v \4 exp[4.0-4.0(arctane)/e] 
0~22- Z 2\ Vth ) l-exp(-27r/e) 

0-21 = 7.42 x 10- 18 (1.66(^)- 205 

v th 

see reference 



■0.66(^)- 



AAZN97 
AAZN97 
AAZN97 
AAZN97 
AAZN97 
AAZN97 
AAZN97 
GJ07 

GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 
GJ07 

AAZN97 
AAZN97 
AAZN97 
AAZN97 
AAZN97 

GJ07 

GJ07 



AAZN97 
AAZN97 
AAZN97 
AAZN97 
AAZN97 
AAZN97 
GJ07 



In the table above, T denotes the gas temperature in Kelvin, T e the gas temperature in eV, T r the temperature of radiation in K. dex(x) = 10 x . 
Acronyms: AAZN97: Abel, Anninos, Zhang, Norman (1997); GJ07: Glover and Jappsen (2007). 
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Table A.2. Heating/cooling rates 



Process 


Rate 


Reference 


Heating/Cooling 


Collisional excitation cooling 


7.5xl(T 19 (l + V75)" 1 exp(-11838/T)n e n H 
9.1xl0" 27 (l + VTsr^^-^expC-llSSS/T^nHe 
5.54 x 10" 17 (1 + V^) _1 T-°- 397 exp(-473638/T)n e n+ e 


AZAN96 


Collisional ionization cooling 


2.18 x 10- n kin e nH 
2.18 x 10~ 11 k 1 n e n H 
3.94 x 10- 11 k 3 n e n H e 
8.72 x 10- n k 5 nen+ e 

5.01 x 10" 27 (1 + V77)- 1 T-° 1687 exp(-55338/T)n 2 n5 e 


AZAN96 


Recombination cooling 


8.7 x 10- 27 T 1/2 T- a2 (l + T^ 7 )- 1 n e n+ 
1.55 x 10- 26 T 03647 n e n+ e 


AZAN96 








i.z4 x iu i [i.u + u.Jexp(— *3>4uuu/ 1 ;jexp(— 4 /uuu/ 1 jn e n He 






j.4o xiu i ^i.u+1^ ) n e n He 

i.4j x iu vi U-i + u.J4exp^— — log^i )) /Jjjn e ^n H + n He +n He j 
5.65 x 10" 36 (1 + z) 4 [T - 2.75(1 + z)]n e 




Bremsstrahlung cooling 


A7A NTQ£» 


Compton cooling 


AZAN96 


Molecular hydrogen cooling 


see reference 


GP98 


HD cooling 


see reference 


LNA05 


H 2 Photodissociation heating 




AAZ97 


Solomon process 


r = 6.4 x 10" 13 H 2 k27 




Direct heating 


r = H 2 JT 4^ 28 (v)^(hv - hv th )dv 

T = 2.93 x 10" 12 k 2 n H + 5.65 x lO^k^n^^^) 




Gas-phase H 2 formation heating 


GJ07 


H 2 Collisional dissociation cooling 


A = 7.2 x 10- 12 (k 9 n H + k 10 n H2 )n H2 


GJ07 


Cooling from electronic states of hydrogen 


see reference 


SSG10 



In the table above, T denotes the gas temperature in Kelvin, T n the gas temperature in 10 n , T r the temperature of radiation in K. dex(x) = 10*. 
Acronyms: AAZN97: Abel, Anninos, Zhang, Norman (1997); GJ07: Glover and Jappsen (2007); GP98: Galli and Palla (1998); AZAN96: 
Anninos, Zhang, Abel, Norman (1996); LNA05: Lipovka, A. et al. (2005); SSG10: Schleicher, Spaans, Glover (2010). 



